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The fluctuation- dissipation relation is usually formulated for a system interacting with a heat bath 
at finite temperature in the context of linear response theory, where only small deviations from the 
mean are considered. We show that for an open quantum system interacting with a non-equilibrium 
environment, where temperature is no longer a valid notion, a fluctuation- dissipation inequality 
exists. Clearly stated, quantum fluctuations are bounded below by quantum dissipation, whereas 
classically the fluctuations can be made to vanish. The lower bound of this inequality is exactly 
satisfied by (zero-temperature) quantum noise and is in accord with the Heisenberg uncertainty 
principle, both in its microscopic origins and its influence upon systems. Moreover, it is shown 
that the non-equilibrium fluctuation-dissipation relation determines the non-equilibrium uncertainty 
relation in the weak-damping limit. 



I. OPEN QUANTUM SYSTEM 



The fluctuation-dissipation relation (FDR) is usually 
formulated for a system interacting with a heat bath 
at finite temperature in the context of linear response 
theory, where only small deviations from the mean are 
considered. We show that for an open quantum system 
interacting with a non-equilibrium environment, where 
temperature is no longer a valid notion, a fluctuation- 
dissipation inequality (FDI) exists. Clearly stated, quan- 
tum fluctuations are bounded below by quantum dissi- 
pation, whereas classically the fluctuations can be made 
to vanish. The lower bound of this inequality is exactly 
satisfied by (zero-temperature) quantum noise and is in 
accord with the Heisenberg uncertainty principle (HUP). 
FDI violating quantum noise can be viewed as arising 
from HUP violating states of the environment and can 
induce HUP violating states in the open system. In 
fact, in the weak-damping limit the non-equilibrium FDR 
(which must satisfy the FDI) precisely determines the 
non-equilibrium uncertainty relation (which must satisfy 
the HUP). 

In the following section we present the necessary back- 
ground material of quantum open systems, wherein we 
formally categorize quantum noise in terms of its time 
dependence, dissipation, and microscopic origin. Read- 
ers familiar with this may skip to our results in the latter 
sections and refer back as needed. In the penultimate 
section we derive the FDI from a microscopic model and 
contrast it to the usual thermal FDR. In the final section 
we work from the other end and motivate the FDI phe- 
nomenologically, but less generally. This result also pro- 
duces the non-equilibrium uncertainty relation for quan- 
tum Brownian motion, which can be contrasted to the 
finite-temperature uncertainty relation [1-4]. 



II. NOISE AND DISSIPATION 

A. Open Systems and Noise 

Consider a quantum system weakly interacting with 
an environment with interaction Hamiltonian: 

Hi = ^ L„ (g) 1„ , (1) 

n 

expanded as a sum of separable operators, where L„ and 
1„ are system and environment operators respectively. 
The environment coupling operators 1„ will typically be 
collective observables of the environment, with depen- 
dence upon very many modes. The system-environment 
interaction will be treated perturbatively and so the cen- 
tral ingredient is the (multivariate) correlation function 
of the environment: 

cw(t,r) = (l„(t)l m (T)) E , (2) 

where l n (t) represents the time-evolving 1„ in the interac- 
tion (Dirac) picture. In the influence functional formal- 
ism [5] for the quantum Brownian model with bilinear 
couplings between the system and its environment [6- 
8] the correlation function appears as the kernel in the 
exponent of a Gaussian influence functional, called the 
influence kernel £ in Refs. [9, 10]. Alternatively, in quan- 
tum state diffusion [11] this kernel takes the explicit role 
of a noise correlation for complex Gaussian noise. The 
influence kernel, or equivalently, the complex correlation 
function, can be written as a sum of two real parts cor- 
responding to the noise and dissipation kernels [9, 10]: 

a(t,r) = v(t,r)+i n{t,r) . (3) 

complex noise noise dissipation 

The noise kernel v appears in the influence kernel as the 
correlation of an ordinary real stochastic source, whereas 
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the dissipation kernel fi alone would produce a purely ho- 
mogeneous (though not positivity preserving in general) 
evolution. These same roles can also be inferred from the 
Heisenberg equations of motion for the system operators 
after integrating the environment dynamics, producing 
the so-called quantum Langevin equation [ ] . Gener- 
alized to nonlinear systems and environments (yet still 
Gaussian in the influence functional), the environment- 
integrated Heisenberg equations of motion can be ex- 
pressed 

s(t) = +i[n, s(t)} + iJ2 M*)> s (*)] U*) > ( 4 ) 

n 

1„ (t) = e„ (*) - 2 T] [ dr ix nm (t, r) L m (r) , (5) 

m - 70 

for the equation of motion of any system operator S(i), 
where £ n (t) is an operator- valued stochastic process with 
symmetrized two-time correlation i/„ m (i,r) and commu- 
tator 2i/j, nm (t, t). The commutator is typically determin- 
istic for quantum Brownian motion. 

Using the notation of Ref. [13] , the second-order master 
equation [14-1G] of the reduced density matrix p can be 
represented in terms of the noise correlation as 

p=[-iH )P }+£ 2 {p}, (6) 

with the second-order contribution given by the opera- 
tion 

C-1{P\ = X! [ L "' P ( A ™™ <> L ™) t _ ( A „m* Lm) p] , 

nm 

(7) 

where the A operators and o product define the second- 
order operators 



(A nm o L m ) (t) — 




(f,r) {G (i,r)L m (r)} , (8) 



given the free system propagator Go(i,r) : p(r) — > p(i). 
One context in which the influence functional, Langevin 
equation, and master equation all work together seam- 
lessly is in the quantum Brownian motion of linear sys- 
tems [17, 18]. In addition to a a quantum Langevin equa- 
tion for non-commuting operators, linearity makes it pos- 
sible in that case to have a Langevin equation for real 
classical stochastic processes from which general quan- 
tum correlation functions and the master equation can 
be exactly derived [IS]. 

B. The Quantum Noise Correlation 

From its microscopic definition, Eq. (2), the environ- 
mental correlation function is Hermitian in the sense of 

a(t,r) = a\r,t) , (9) 

and also positive definite in the sense of 

/ d n f d7*ft(Ti)a(Ti,7a)f(7s) >0, (10) 
Jq Jo 



for all vector functions f (t) indexed by the noise. Positiv- 
ity and the noise decomposition (3) are the key properties 
from which the FDI arises. In next section we proceed to 
categorize the relevant time-dependence and correspond- 
ing microscopic origins of the environmental correlations. 
Then we will relate these features to the categorization 
of environments into resistive, amplifying, and indefinite. 

C. Categorization of Noise 

Stationary correlations are defined by their invariance 
under time translations, 

a(t,T) = a(t-T), (11) 

and can produce asymptotically stationary (time- 
constant) master equations. Such correlations are pro- 
duced when the environment is in an initially stationary 
state and its coupling operators in the Schrodinger pic- 
ture are constant in time. 

P E (o) = J2M^)hM, (12) 

i 

yielding the correlation function 

a„ ro (t, t) = ^p E (£i) (£i\ 1„ \e i} ) lm |e lJ -) e +l£ ^ T) , 

(13) 

where = £j — £j, \e) denotes the energy basis of the 
environment and pe(s) are its stationary probabilities at 
the initial time. The accompanying characteristic func- 
tion can be obtained quite directly from the mode sum 

a(t) = —J_ duje + ^&(uj), (14) 

yielding the Fourier transform 
a nm (u>) oc27T^p E (£ 4 ) fa] In \£i ~ w) (e,| l m |e, - u) , 

where the underscored proportionality here is strictly in 
reference to the continuum limit of the reservoir which 
relates environmental mode sums to integrals given the 
infinitesimal strength of individual environmental mode 
couplings. This can be more rigorously defined through 
the use of a finite spectral density function in place of 
the infinitesimal environment couplings. 

Also of note are quasi- stationary correlations of the 
form 

a(t,r) = a s (t - r) + Sa(t + r) , (16) 

where as(t— r) denotes a stationary correlation function, 
or, more specifically, the stationary projection of a(t, r), 
while Sa(t-\-r) is an additional non-stationary contribu- 
tion. Such correlations will result from (constant) linear 
coupling to an environment with non-stationary initial 
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state, such as a squeezed thermal state. In these cases 
the stationary projection of the correlation function does 
correspond to the stationary projection of the initial state 
of the environment. As for the non-stationary contri- 
butions, due to their highly oscillatory behavior in the 
late-time limit they typically lose effect asymptotically. 
Therefore, quasi-stationary correlations can produce an 
asymptotically stationary master equation with equiv- 
alent asymptotics as generated by their corresponding 
stationary correlation. 

Finally one can also consider cy do -stationary and 
quasi- cy do -stationary correlations, which exhibit a peri- 
odic time-translation invariance and can produce asymp- 
totically cyclo-stationary master equations. Much of this 
work can be easily generalized to these cases, with the ex- 
ception of detailed balance which cannot be maintained. 



D. Noise Decomposition 

Motivated by the influence functional formalism and 
Heisenberg equations of motion, correlation functions of 
second order can always be decomposed into a real noise 
kernel and dissipation kernel as in Eq. (3): The Hcrmitic- 
ity stated in Eq. (9) leads to the relations 
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u(t,r) = - [a{t,T) + a T (T,t)] , 
H(t,r) = i [a(t,r) - a T (T,t)] . 



(17) 
(18) 



The two kernels naturally decompose the second-order 
operators of Eq. (8) into their Hermitian and anti- 
Hermitian parts, in the ordinary sense of Hilbert space 
operators: 



N, 



-l Mr. 



(19) 



diffusion dissipation 

(N„ m ©L m ) = f dTis nm (t,T) {G (<,r)L m (r)} ,(20) 
Jo 

(M nm oL m ) = / dT[i mn (t,T) {Go(t,r)L m (r)} (21) 
Jo 

The second-order master equation Ci{p\ can then be 
expressed entirely in terms of Hermitian operators as 

- ^ [L„, i{(M mn o L,„), p} + [(N nm o L m ), p]] . (22) 

nm 

Here the noise coefficients describe diffusion while the 
dissipation coefficients describe dissipation (or amplifica- 
tion), renormalization and other homogeneous dynamics. 

The correlation function a(t, t) is positive definite and, 
therefore, the noise kernel v>{t 1 r) must also be positive 
definite. The dissipation kernel fi(t, r) is not positive 
definite, but it is related to the damping kernel 7(i, r), 
which is given by 



d_ 



and can be positive definite, negative definite, or indefi- 
nite. The dissipation kernel coefficients can then be ex- 
pressed in terms of damping kernel coefficients 

(M nro o L m ) = (r„ m o L m )(i) - 7„ m (i, t) L m (t) 

s, v ✓ v v s ^ v ^ 

dissipation damping renormalization 

+ 7nm (i,0){G o (i)L ro (G)} (24) 

S * ' 

slip 

(r„ m oL m ) = J dTj nm {t,T) |G (i,r)L m (r)| (25) 

L m (t) = +i [H(i) , L m (t)] + ^L ra (t) , (26) 

which can be used to place Eq. (22) into a form much 
like the QBM master equation [7, 18, 19]. The slip is a 
transient effect as a result of dealing with a factorized 
initial state, which is modified if one considers correlated 
initial states. The renormalization is a permanent shift 
of the system Hamiltonian which would diverge in the 
limit of local damping and is typically canceled with a 
counterterm introduced with the bare system action or 
the system-environment interaction. 

For stationary correlations a(t — t) with characteristic 
function (Fourier transform) a(uj) = J 00 dr e~ lUT 'a(r) , 
the noise and damping kernels are then Hermitian in both 
noise index and frequency argument 



7M = 7 f M = 7*(~w) . 
v{w) = if '(u) = £>*(— u) , 



(27) 
(28) 



and by Bochner's theorem both a(uj) and u{<jS) are 
positive-definite for all frequencies. Again the damping 
kernel 7(w) may be positive definite, negative definite, or 
indefinite. 



E. Classification of Damping Kernels 

Environments with positive-definite damping kernels 
are damping or resistive environments, while those with 
negative-definite damping kernels are amplifying. If the 
system coupling operators L„ are position operators, the 
damping terms correspond to forces linear in momentum. 
Stationary correlations are the easiest to dissect and the 
most well behaved. Their dissipation and damping ker- 
nels are related by 



(29) 



and from the definition of the dissipation kernel in 
Eq. (18) and the double Hermiticity in Eq. (27)-(28), the 
damping kernel will be most-generally positive or nega- 
tive definite if we have a strict inequality between positive 
and negative energy argumented environmental correla- 
tions. 



7(*,t), 



(23) 



d(-M)>a*(+H) 
a(-M)<a*(+M) 



(Damping) , 
(Amplifying) . 



(30) 
(31) 
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From Eq. (15), one can show that damping environments 
result when the initial stationary probability of the en- 
vironment Pe(e) is a monotonically decreasing function 
of the environment energy. Amplifying environments re- 
sult from monotonically increasing functions or popula- 
tion inversion. The most common example of each being 
positive and negative temperature reservoirs. 

Given our damping representation of the master equa- 
tion coefficients, one can determine the dynamics of the 
system energy from the super-adjoint of the master equa- 
tion [201: 



Ln, (N B ra O L m ) 



+ ■ 



(32) 

where we have neglected any power generated by the slip 
and time-dependence intrinsic to the coupling. Using 
the zeroth-order solution p(t) = G (t) p(0) and symme- 
tries of the damping kernel, the second-order expectation 
value for the cumulative energy lost through damping is 
given by 



dr 2 7„ 



(ri,r 2 )Tr L n (n)p(0)L m (T 2 ) 



which will be strictly dissipative for a positive-definite 
damping kernel. This expression also contrasts nonlo- 
cal damping to local damping. Delta correlated damping 
kernels are strictly dissipative at every instant of time 
whereas nonlocal damping kernels are only assured to be 
net dissipative in the full time accumulation. 



III. NON-EQUILIBRIUM RELATIONS 

A. Non- Equilibrium Fluctuation-Dissipation 
Relation & Inequality 

From the definitions of the multivariate noise kernel 
u{ijS) and damping kernel 7(w), Eqs. (3)-(18), one can 
prove the fluctuation-dissipation inequality: 



v{uS) > ± u) 7(w) , 



(34) 



here for stationary correlations, and in the Fourier do- 
main where the u would denote energy level transitions 
of the system. To prove this one simply notes that the 
noise kernel is the sum of two positive-dehnite kernels 
whereas the dissipation kernel is given by their differ- 
ence. The essential point is that if there is any damp- 
ing, or amplification, there will be quantum noise and 
Eq. (34) determines its lower bound. This is quite a de- 
parture from classical physics where noise can be made to 
vanish in the zero-temperature limit, although the lower 
bound of this inequality is satisfied by zero-temperature 
quantum noise since ct(|a;|) = in that case. 

For the case of one collective system coupling, coupled 
to one or more environments, it is sufficient to define a 



fluctuation-dissipation relation 



v(u) = k(uj) 7 (w) , k(u) = V !j J \ , 



(35) 



with k(w) the fluctuation-dissipation kernel [7, 8] which 
relates fluctuations to dissipation. For multivariate noise 
one might use the symmetrized product 



£>(w) 



1 



\k(u) 7(w) + 7(0;) k(uj)] , 



(36) 



which would ensure k(uj) to be positive definite if 7(w) 
is, in accord with this being a (continuous) Lyapunov 
equation [21]. We will use this particular definition for 
quantum Brownian motion in the next section. Inequal- 
ity (34) can now be restated as 



k(ui) > \ui\ 



(37) 



for damping environments. Typically k(uj) will contain 
dependence upon the precise nature of the environment 
couplings 1„. 



B. Equilibrium Fluctuation-Dissipation Relation 

Let us consider a constant system-environment interac- 
tion, constant environment Hamiltonian, and initial sta- 
tionary probabilities of the environment given by pe(^)- 
If the FDR is to be independent of precisely how the sys- 
tem and environment are coupled, then one can work out 
from the microscopic theory that the FDR kernel must 
be a scalar quantity, directly related to the initial state 
of the environment by way of 



«M _ Pe(£ - u) +Pe(e) 
ui Pe(e - oj) -Pe(e) 



(38) 



for all e. To prove this one first applies relation (15) to 
definitions (17)-(18), and notes that if the dissipation and 
noise are related in a manner independent of the coupling 
then the two kernels must be related term-by-term in a 
sum over couplings. 

But such an equality between the FDR kernel and 
mode probabilities implies the functional relation 



1 



1 



(39) 



where the lj translations can factor out. This factoriza- 
tion property is unique to exponential functions; there- 
fore, only the thermal distribution pe(e) oc e~^ £ can pro- 
duce a fluctuation-dissipation relation which is generally 
coupling independent. We then have that 



kt{lo) = lu coth^ 



27V 



(40) 



for the thermal distribution. One should be careful to 
note that the thermal FDR is not special because it ex- 
ists, nor because of its simple form, but because of its 
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invariance to couplings (model-independence) . In a more 
general context, the thermal FDR is also special because 
it ensures a relaxation to detailed balance in a coupling- 
invariant manner. In fact, these properties can be shown 
to be equivalent [13]. 

As a concrete example of an elegant yet non- 
equilibrium FDR, the late-time dominating stationary 
correlations for linear coupling to a squeezed thermal 
reservoir [22, 23] will produce the FDR kernel 



= cosh[2r(w)] wcoth^ 



—) 
277 



(41) 



where r(to) is the squeezing parameter, which may be 
allowed to vary with the energy scale. One can easily see 
that this FDR also obeys inequality (34) as it must. 



IV. NON-EQUILIBRIUM UNCERTAINTY 
PRINCIPLE 



In the context of second-order perturbation theory, 
quantum noise is effectively Gaussian in the influence 
functional, and Gaussian noise is equivalent to that aris- 
ing from linear coupling to a bath of harmonic oscillators. 
Therefore any violations of Eq. (34) must correspond to 
environment oscillators in a non-quantum state. In the 
phase-space or Wigner function representation [ ] , HUP 
violating states of the environment can be constructed 
which violate the quantum FDI. Such is the case for 
the classical vacuum, which has vanishing fluctuations 
yet finite damping. Now we shall show that FDI violat- 
ing noise can also relax the system into a HUP violating 
state. 

Let us consider weakly influencing a system of oscilla- 
tors at resonance, all with mass m and frequency w, via 
position-position coupling to some phenomenological set 
of noise processes with resistive correlation a(uj). We do 
not assume the system-environment couplings to be iden- 
tical, nor do we neglect the presence of cross-correlations 
among the noise processes. 

From the results of Ref. [18, 19] and the second-order 
master equation coefficients (7), the damping kernel 7(w) 
will play the role of the dissipation coefficients and the 
noise kernel 0(uj) will play the role of the normal dif- 
fusion coefficients in the Fokker-Planck or master equa- 
tion. Integrals over the two kernels will then provide 
the system renormalization and anti-diffusion coefficients 
respectively. Given sufficient dissipation and bandwidth- 
limited correlations, the system will relax into a Gaussian 
state which satisfies the Lyapunov equation 



(42) 



for the momentum covariance, which has elements (piPj), 
and 

Vxx = 7 ^ a pp > (43) 

(T XP = 0, (44) 



for the remaining covariances in the phase-space (Wigner 
function) representation [ ], and to lowest order in the 
system-environment interaction. Comparing Eq. (42) to 
Eq. (36), we can express our covariances 



1 



'pp 



2muj 
m 



(45) 
(46) 



in terms of the FDR kernel k(u>). So far our FDR ker- 
nel remains phenomenological and not microscopically 
derived, however it must be positive definite for this to 
describe a physical state. As k(lj) is positive definite, we 
may transform to the basis in which it is diagonalized. If 
k(lo) is a scalar quantity, then this is simply the normal 
basis of the free system. For each mode in this basis we 
have the phase-space covariance 



2mu 



K n {U) 









(47) 



which must then satisfy the generalized Heisenberg un- 
certainty relation due to Schrodinger [25, 26]: 

(Ax 2 )(Ap 2 )-±({Ax,Ap})>\, (48) 
or in terms of the phase-space covariance 

det(o-) > - (49) 
and, therefore, it must be the case that 

> w , (50) 

for all u> > 0. But this is equivalent to our previous 
statement 



k(ui) > lu , 



(51) 



in terms of positive definiteness as ui is a scalar quantity. 
So not only do FDI violating correlations arise from HUP 
violating states, they can also produce HUP violating 
states via dissipation and diffusion (and decoherence). 

Furthermore we can say that in the weak-damping 
limit, the scalar FDR kernel k(ui) precisely determines 
the (asymptotic) non-equilibrium uncertainty product 



det(er) 



1 k(cu) 

2 LU 



(52) 



for a single system mode of frequency lo. Larger FDR 
kernels naturally produce larger uncertainty and mini- 
mal FDR kernels (zero temperature) produce minimal 
uncertainty. Non-perturbative results require evaluation 
of the exact expressions found in Refs. [18, 19] for a single 
system oscillator and multiple system oscillators respec- 
tively. 
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V. DISCUSSION 



In this paper we have derived a fluctuation-dissipation 
inequality (FDI) for an open quantum system inter- 
acting with a non-equilibrium environment from the 
microscopically-derived environment correlation function 
and recovered the well-known fluctuation-dissipation re- 
lation (FDR) for a thermal environment. The FDI is 
a very general statement contrasting quantum noise to 
classical noise, and is satisfied even for non-equilibrium 
fluctuations. Simply put, quantum fluctuations are 
bounded below by quantum dissipation, whereas clas- 
sically the fluctuations can be made to vanish. The 
lower bound of this inequality is exactly satisfied by 
zero-temperature noise and is in accord with the Heisen- 



berg uncertainty principle (HUP). FDI violating corre- 
lations arise from HUP violating states of the environ- 
ment and can relax the open system into HUP violating 
states. Therefore, the FDI can be viewed as an open- 
system corollary to the HUP both from microscopic and 
phenomenological considerations. Analogously, the non- 
equilibrium FDR also determines the non-equilibrium 
uncertainty product, most directly in the limit of weak- 
damping [cf. Eq. (52)], and the corresponding FDI im- 
plies the HUP. 
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